Spin dependent non-tunneling transport in granular ferromagnet 
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We have shown that, at low temperature, the effect of localization of charge carriers in granular 
metals cannot be ignored. The electron localization in granular ferromagnetic metals can give rise 
to non-linear transport in the moderate electric field regime, which is dominated by spin dependent 
inter-grain non-tunneling emission. In contrast to the general opinion that spin dependent tunneling 
gives rise to low field magnetoresistance in granular ferromagnetic metals we observe that spin 
dependent non-tunneling transport is primarily responsible for the sharp magnetic field dependence 
of resistance, at least in the low temperature and moderate electric field regime. 

PACS numbers: 72.25.-b, 73.63.-b, 85.75.-d 



Recent emergence of spin based electronics or spin- 
tronics has initiated extensive studies on spin polarized 
transport in magnetic tunnel junctions, granular ferro- 
magnetic films and nanocomposites. The low field mag- 
netoresistance in magnetic tunnel junctions (MTJ), and 
granular ferromagnetic systems is attributed to spin de- 
pendent tunneling 0.0. 0]. However, it has been shown 
recently that the superposition of tunneling conduction 
and ballistic channels can lead to inversion of magne- 
toresistance 4t] in MTJs. This letter explores the possi- 
bility of whether there exists any transport mechanism 
other than spin dependent tunneling which might af- 
fect low field magnetoresistance in granular ferromag- 
nets. Although the study of granular ferromagnetic met- 
als spans more than three decades Q, recently, there is 
a renewed interest in this field with the granular sys- 
tems based on "half metallic" materials like CrC>2 and 
Lao.67SrnQ3Mn03 showing enhanced low field magnetore- 
sistance [ijj leading to the possibility of their application 
as magnetic field sensors. We will analyze the effect of 
electron localization in granular ferromagnetic systems at 
sufficiently low temperature. The consequent dominance 
of spin dependent inter-grain non-tunneling emission over 
inter-grain tunneling in the moderate electric field regime 
will be discussed. 

For the preparation of granular film with almost uni- 
form grain size distribution and quasi-homogeneous dis- 
tribution of conducting and non-conducting regions, we 
have utilized the well studied strain relaxation properties 
of manganite films. Manganites like Lao.67Sro.33MnC>3 
(LSMO) can be grown epitaxially on single crystal in- 
sulating substrate LaAlC>3 (LAO). However, due to fi- 
nite lattice mismatch ({cilsaio — &lao) I &lao ~ +2%) 
with the substrate, thin films of Lao.67Sro.33Mn03 un- 
dergoes a biaxial in plane compressive stress leading to 
higher out ofplane lattice parameter compared to the 
substrate |5j, 111 III ■ Relatively thick films (usually above 
50 nm) become strain relaxed forming granular nanos- 
tructures and the grain size depends on the thickness of 
the film. In granular manganite systems, it is an estab- 
lished fact that the core of grain is ferro-magnetically 
ordered but the grain surface is spin-disordered and in- 



sulating [1,0. 

Thick films of LSMO, having thickness ~ 3000A, were 
grown on LaAlO3(100), using pulsed laser deposition. 
The surface morphology of the film was inspected, using 
atomic force microscopy which shows orderly arrange- 
ment of connected grains with average size 80 nm (Fig:^|. 
In order to prove that this is not a mere surface corru- 
gation, the x-ray diffraction 9 — 29 scan around (002) 
peaks of LSMO/LAO are shown in Fig:|2K and the mag- 
netic measurements in Fig: The Bragg peak corre- 
sponding to that portion of the film which matches with 
the in-plane lattice parameters of the substrate leading 
to higher c-axis parameter is indicated in the figure as 
(1) which gradually relaxes to the bulk lattice param- 
eter value shown by the higher angle peak among the 
twin peaks separated by a plateau indicated as (2). The 
intensity of Bragg peak which corresponds to strain re- 
laxed part of the film is much higher compared to the 
strained part indicating that the granularity extends deep 
inside the film and hence bound to influence the transport 
and magnetic properties of the film. The magnetic char- 
acterization was done using a Quantum Design SQUID 
magnetometer which shows the ferromagnetic transition 
temperature T c at around 360 K. The large irreversibil- 
ity in the temperature dependence of zero field cooled 
(ZFC) and field cooled (FC) magnetization (M) together 
with the hysteric M — H loop, where the magnetization 
does not saturate even at 2 kOe, are typical of granu- 
lar ferromagnets (Fig: |2K,B). The virgin M — H curve 
shows (Fig:|2J) an 'S'-like shape with the point of inflexion 
around 100 Oe indicating non-ferromagnetic inter-grain 
coupling. The coercivity of the film has roughly the same 
value. 

The transport properties were studied using four probe 
technique with the magnetic field being in-plane to the 
electric field. The granular film exhibits a sharp rise in 
resistance (Fig:0j3) with lowering of temperature below 
around 6 K in the absence of external magnetic field. At 
higher bias current, the rise in resistance reduces drasti- 
cally suggesting non-ohmic transport at low temperature. 
A significant reduction of percentage rise in resistance on 
application of magnetic field is observed. Hence Ziese's 
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FIG. 1: AFM picture of the strain relaxed film over 2 x 2/wn 2 
area showing the granular surface morphology along with a 
schematic cross-sectional blown-up view with strain-relaxed 
granular part and a thin strained bottom part, the grains 
(the ferromagnetic ally ordered core part) being assimilated to 
localized sites at low temperature, leading to the dominance 
of inter-grain non-tunneling transport. 
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FIG. 2: A: X-ray diffraction pattern of the LSMO film on 
LAO substrate. The c-axis parameter of the film relaxes 
gradually to the bulk value of LSMO. The broad hump at 
the lower angle (1) indicates the strained part whereas the 
two broad peaks (2) and the plateau in between corresponds 
to the gradual relaxation of the strain. B: The conduction is 
due to the combined effect of tunneling through the barrier 
and thermal activation over the barrier. The effective barrier 
height depends on the relative spin orientation of adjacent 
sites. The conduction is in the moderate electric field regime 
where R, - R^ 1 « R, - R, 



model of quantum interference effect ^(j is not applica- 
ble here. Below 6 K, the strain relaxed film exhibits low 
field magnetoresistance which is extremely sensitive to 
the applied bias (Fig:0p). 

The electric field dependence of conductivity at 2 K, 
plotted in Fig: 0)3, shows that the logarithm of conductiv- 
ity is linearly proportional to square root of the electric 
field (E). At 5 K, very close to the point of upturn in re- 
sistance with lowering of temperature, the transport be- 
comes nearly ohmic (Fig:0p). Such electric field depen- 
dence is strikingly similar to the charge tr ansp ort mecha- 
nism called Poole-Frenkel emission (PFE) 11]. The PFE, 
usually observed in disordered semiconducting systems, 
is a non-tunneling transport phenomenon where the 
charge carriers are thermally activated over the poten- 
tial barrier from localized charge traps in the presence of 
electric field. The reduction in barrier height cjypp is pro- 
portional to the square root of applied electric field and 
the conductivity is given by G = Go exp (PppE 1 / 2 /KT) , 



where (3pf is a constant. 

However, in this case what is fundamentally different 
from traditionally observed PFE is the tunability of the 
phenomenon by application of even a small amount of 
magnetic field. It is evident from Fig:|3J\ that the thresh- 
old electric field for the onset of charge emission E t h (in- 
dicated by the vertical arrow) and the slope of the curve 
are higher at 100 Oe magnetic field compared to that 
at 1 kOe. The origin of PF type emission in granular 
metal can be explained following Mostefa's model [l2| . 
This model considers a granular metallic system with a 
narrow grain size distribution with almost equal Fermi 
levels for all grains and as a result the activation energy 
is entirely due to the potential barrier between two ad- 
jacent grains. In each grain, electrons can be considered 
as being localized at sufficiently low temperature. At 
low temperature, the leakage of electron wave function 
outside the grain is minimal, which is particularly true 
for insulating grain surfaces, as in this case This 
justifies the assumption of grains as localized sites at suf- 
ficiently low temperature. In the weak field regime, the 
transport is purely thermally activated. The electric field 
will affect two important parameters of the system if it 
is above the weak field regime. 

Firstly, the shape of the potential barriers and hence the 
decay of the wave function associated with the electronic 
states. The tunneling probability P(Rij) between two 
grains assumed to be located at sites Ri and Rj separated 
by a distance Rij = \Rj — Ri\ can be written, using WKB 
approximation, 



P(Rij) -exp 



{2m[V{x) -e l ]} 1/2 dx 



Ri 



Here V{x) is barrier height, x is the position of the elec- 
tron on the path R^ and Si is the electron energy in the 
grain located at ith site. It is assumed that the potential 
barrier is in the moderate field regime where the electron 
is emitted as a free particle into the insulator over a path 
much shorter than the total distance between the grains. 
The effective barrier height is 4>i(x) — V(x) — E; L . Taking 
these factors into consideration and assuming directed 
flow of electrons we get, P(Rij) — exp { — where A 
is a constant. 

Secondly, the activation energy E a required for one 
electron to hop from a neutral grain to another neutral 
grain. In the non-negligible field regime, the activation 
energy required to transfer an electron from site i to site 
j is, E a rj — Sij — eE.Rij , where, eij = Ej — Ej . Hence the 
activation probability, 



KT 



eERij cos 8 
KT 



The electrical conductance G%j between sites i and j 
can be calculated from the total transport probability Pij , 
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FIG. 3: A: The zero field cooled (ZFC) (open symbols) and 
field cooled (FC) (filled symbols) M vs. T curves at magnetic 
fields H = 20 Oe and H — 100 Oe showing irreversibility 
between the ZFC and FC curves which is an evidence for the 
granularity of the film. B: Virgin M vs. H curve for the strain 
relaxed film at 2 K, exhibiting inflexion at 100 Oe shown by 
the slope of the curve plotted as continuous line. C: M vs. H 
hysteresis loops at 2 K, showing coercivity of about 100 Oe. 
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FIG. 5: Sensitivity of the PF type emission to magnetic field. 
Plot of logarithm of dc conductivity vs. E 1 ^ 2 in the antiparal- 
lel and parallel magnetic configuration with the corresponding 
magnetic field dependence of magnetoresistance at the same 
temperature in the inset for the strain relaxed film (A) and 
(LSMO)i_ a: (ALO) a: composite with x = 0.5 (B). The thresh- 
old electric field E t h is indicated by the vertical arrow. The 
additional inset (II) in (B) describes the temperature depen- 
dence of resistance of the same sample with a rapid increase 
in slope below 10 K shown by the continuous line. 



FIG. 4: A: Rise in resistance with lowering of temperature 
below 6 K at different bias currents, suggesting non-ohmic 
transport below 6 K. Inset: Magnetic field dependence of MR 
at 3 K at different bias currents. B: Plot of logarithm of dc 
conductivity vs. E 1 ^ 2 in absence of magnetic field at 2 K 
and 5 K for the same sample. The linear region in the curve 
(for 2 K) indicates thermally assisted non-tunneling emission 
due to localization of charge carriers at low temperature. The 
non-tunneling emission is almost suppressed at 5 K. 

which is given by the product of the tunneling probability 
across the barrier and the activation probability over the 
barrier. Assuming that the temperature is low enough 
to help the hopping transport, the conductance can be 
calculated employing percolative treatment. 

However, Mostefa's model alone cannot explain 
the sharp magnetic field dependence discussed earlier. 
Therefore, we incorporate an additional inter-grain mag- 
netic exchange energy (E m ) term in the model. It arises 
due to the relative spin orientation of the two adjacent lo- 
calized sites when the spin is conserved in hopping from 



one site to another. The value of E m depends on the 
overlap of electronic wave function across the grains. If 
the spin orientation is anti-parallel an additional energy 
+E m is necessary for charge carrier emission while less 
energy — E m is required if the orientation is parallel. Ob- 
viously, the value of E m can easily be manipulated by 
applying even a tiny magnetic field H . Hence the pre- 
scription is to simply replace the effective potential bar- 
rier 4>i{x) by 4*™{x) = 4>i(x) ±Ell(H) and the activation 
energy E a ij by E a ij ' m = e 4j - eE.R^ ± E»(H). Now we 
employ the percolation criterion Gy > Gc, where Gc 
is the critical conductance for percolation. This sets a 
limit on the energy and the length scales (sij < Em and 
Rn < Rm) involved and one arrives at, III _ -5m. < i Fa 

where, e M = (if^) Y\ S M = (^) Y' and Y' = In ^ 
Here the magnetic exchange energy has been expressed in 
the unit of as E m = ^ m £ij, where £ m is a dimension- 
less quantity. Then following Mostefa's derivation one 
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obtains, 




where -BfX TT ^ VI i 6m f° r a given system, with || 
and 1"f indicating the antiparallel and parallel magnetic 
configuration between nearest neighboring grains, respec- 
tively. Evidently, the value of B is higher in the antipar- 
allel configuration compared to the parallel one. 

It is interesting to compare the magnetoresistance (In- 
set, Fig: EK) with the conductivity curves at different 
magnetic fields in Fig: at low temperature. The mag- 
netoresistance is positive at 100 Oe, which is equal to the 
coercivity and the inflexion point in the virgin M — H 
curve (Fig: OJ). This suggests that at this point, the 
spin orientation between the adjacent grains becomes 
anti-parallel, resulting in increase of exchange energy. 
The increase in the slope of the logarithmic conductiv- 
ity vs. E 1 ' 2 curve at H = 100 Oe in Fig: [5] compared 
to that at if = 1 kOe is consistent with our model. 
The value of exchange energy can be estimated from 
the slope of the curves at 100 Oe (antiparallel config- 
uration) and 1 kOe (parallel configuration) using the re- 
lation £ m = (B n 2 - B n 2 )/(B n 2 + B TT 2 ) which gives 
£ OT = 0.79 and therefore E m = 0.79£y. The validity 
of the non-tunneling emission being limited to the mod- 
erate electric field regime (as defined in Fig:|23) is deter- 
mined partly by the potential barrier height. Therefore, 
the shift in the threshold electric field (Eth) for the onset 
of PF type emission should be magnetic field dependent 
since magnetic exchange energy contributes substantially 
to the total barrier height. In the antiparallel configura- 
tion, the value of E t h should be higher in comparison 
with the parallel configuration, which is consistent with 
our observation. A rough estimation of the average total 
activation energy (including the magnetic contribution) 
is possible from the fitting of temperature dependence 
of resistance (at low temperature and bias) in absence 
of magnetic field by an activation term with a metallic 
contribution in parallel. The metallic channel includes 
whatever the effect of the thin continuous strained part 
might have. Consequently, the magnetic exchange en- 
ergy turns out to be E m ~ 0.24 K. Similar value (~ 0.23 
K) can be extracted from the difference in the total ac- 
tivation energy values in the absence and in presence of 
magnetic field, respectively. 

In order to find out whether the same phenomenon is 
operative in a different micro-structure, we have stud- 
ied magnetotransport properties in a nanocrystalline 
LSMO/ A^C^ALO) granular composite with different 
weight concentrations x of ^203 (where x — — 0.5). 
The average grain size of the sol-gel prepared LSMO 
was measured, using Transmission Electron Microscope 
(TEM) and Scanning electron microscope (SEM), which 
turns out to be about 50 nm. Fig: |3j3 shows similar 
spin dependent non-tunneling emission in LSMO/ALO 



composite with x = 0.5 at 3 K. The temperature de- 
pendence of resistance shows a low temperature upturn 
with a rapid increase in the slope below 10 K (Inset II, 
Fig: 03). All other samples with different values of x, 
exhibit similar behavior. This suggests that the effect 
of localization at low temperature and the consequent 
dominance of non-tunneling effect might be generic to 
a large variety of granular ferromagnetic metals. Em- 
ploying the same method for estimation of magnetic ex- 
change energy described in the previous paragraph we 
obtain £ m = 0.04 only. The temperature dependence of 
resistance can be fitted using the parallel channel con- 
duction model of ref. 0] , which gives a rough estimate 
of activation energy to be ~ 0.8 K. The reduction in the 
slope of the curve as well as E t h with applied magnetic 
field is not as pronounced as in the case of strain-relaxed 
film indicating that the value of exchange energy should 
be comparatively smaller. The small value of exchange 
energy in this case can be understood considering the 
larger average inter-particle separation in the composite. 
Thus the sensitivity of the slope of the logarithmic con- 
ductance vs. E 1 / 2 curves in the non-tunneling regime 
as well as that of E t h to the magnetic field depends on 
the relative contribution of the magnetic exchange to the 
total activation energy. This factor is determined by the 
microstructurc of the granular system along with the spin 
polarization of the ferromagnetic grains. 

To conclude, we have observed spin dependent 
non-tunneling transport in LSMO granular film and 
LSMO/ALO metal-insulator nanocomposites, a phe- 
nomenon which appears to be generic to ferromagnetic 
granular metallic systems, particularly at low tempera- 
ture. In contrast to the usual Poole-Frenkel emission, 
the electric field dependence of conductivity is sensitive 
to external magnetic field. The observed phenomenon 
has been analyzed considering the dominance of ther- 
mally assisted inter-grain non-tunneling transport (with 
the grains acting as localized sites) in the moderate elec- 
tric field regime and the existence of an additional inter- 
grain magnetic exchange barrier. 
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